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BROWNIAN MOTION. 



NATHANAEL ENRIQUEZ 

Abstract. In this work we introduce correlated random walks on Z. When 
picking suitably at random the coefficient of correlation, and taking the average 
over a large number of walks, we obtain a discrete Gaussian process, whose 
scaling limit is the fractional Brownian motion. We have to use two radically 
different models for both cases ^ < H < 1 and < H < |. This result 
provides an algorithm for the simulation of the fractional Brownian motion, 
which appears to be quite efficient. 



1. Introduction 

The fractional Brownian motion appears to be a very natural object, for its three 
fundamental and characteristic features of being a continuous Gaussian process, 
self-similar and with stationary increments. By self-similar, we mean that there 
exists a real number H e]0,l[, such that the finite-dimensional distributions of 
{T-"BH{Tt), i > 0} do not depend on T. The parameter H is called the Hurst 
parameter or the index of self-similarity. Under these conditions, it is not hard to 
check that the covariance function must have the following form: 

E[BHis)BHit)] - ^{s'" + t'" ~\s~ t\-'"). 

A major consequence of this fact is that if we introduce the sequence of increments 
GhU) = Bnij) — BhU — 1), j — 1, 2, also called fractional Gaussian noise, we 
note that they are strongly correlated (for H ^ More precisely, 

E[GH{j)GH{j+k)] - H{2H-l)k^"-^ (1) 

k — *oo 

We denote the two radically different behaviours on both sides of ^: for H < 1/2 the 
increments are all negatively correlated, which corresponds to a chaotic behavior, 
whereas for H > 1/2 the positive correlation between the increments corresponds to 
a more disciplined behaviour. We refer to Chapter 7 of Samorodnitsky and Taqqu 
for more information about fractional Brownian motion. 
There is a well-known representation introduced by Mandelbrot and Van Ness 
m of the fractional Brownian motion as an integral of a kernel function with re- 
spect to the usual Brownian motion. The method of approximation that uses the 
discretization of this integral leads to an algorithm obliging to store a lot of data 
in memory, and to deal with non smooth functions in the case Q < H < ^. This 
inconvenience is nicely discussed by Carmona and Coutin in Q , and they propose 
a way to reduce it. 
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The construction we propose is based on correlated random walks: it consists 
in discrete processes such that the law of each move is ruled by the value of the 
previous move. We refer to Q and the references therein for more information 
about these processes. We note that the decay of correlation for such processes is 
exponential, but that a mixture of these processes, adapted to the value of the index 
of self-similarity, leads to walks whose correlations satisfy (1). The superposition of 
a large number of such walks yields a discrete Gaussian process whose correlations 
fulfill the conditions of Taqqu , so that its scaling limit is the fractional Brownian 
motion. We have to distinguish both cases ^ < H < 1 and < H < ^. In the 
first case, the equivalence relation (1) is the only condition to check. The case 
< -ff < ^ is more delicate: a compensation relation between the correlations 
has to be satisfied simultanously. That leads us to introduce two different types of 
correlated random walks for both cases. 

This construction reminds the constructions using renewal processes of Mandel- 
brot § and Taqqu and Levy It reminds even more the construction coming 
from traffic modeling by Taqqu, Willinger and Sherman |p^ . To our knowledge, all 
these methods are restricted to the case ^ < H < 1. In that case, we shall discuss 
the differences between our construction and 1 14 . Finally, we want to mention that 
random media has already been used by Kesten and Spitzer |^ to get convergences 
towards certain self-similar processes with index of self-similarity bigger than i . 

In the last part, we discuss the speed of the convergence stated in the previous 
parts. We indicate how to deduce a quite fast algorithm that does not ask to keep 
a lot of datas in memory. 



2. The case 1/2 < iJ < 1 

2.1. The correlated random walk. We introduce first our basic tool: the cor- 
related random walk with persistence p. It is a process evolving on Z by jumps of 
-|-1 or -1, whose probability of making the same jump as the previous one is p: 

Definition 1. For any p e [0, 1], the correlated random walk with persistence 
p is a "L-valued discrete process, such that : 

_ Xl = 0, P{Xl = -1) = 1/2, P{Xl = 1) = 1/2. 

_ Vn > 1, eP := — X^_-^ equals 1 or -1 a.s. 

- Vn > 1, P(4+i = el\<7{Xl 0<k<n))^p. 

This process is not Markovian, but if we define a state as the position of the 
process on Z, coupled with the sign of its last jump, we have to deal with a Markov 
process onZx{— 1,1}. In fact, this process consists in alternate falls and rises of 
i.i.d. geometric length with parameter p. 

We can compute the correlations between two steps distant from n: 

Proposition 1. Vm > l,n > 0, £^[e^e^,+„] = {2p - 1)" 
Proof: Vn > 1, 

E[eP+,l,.^,\a{Xk,0< k < n)] + E[ei+,l,P^_MXk,0 < k < n)] = 
(2p- l)l,P^l - {2p- l)le^„=-i = (2p- 1)£P. 
Consequently, conditioning by cr(Xfe, < k < m + n), we get 
Vm > l,n > 0, i?[C+„+i£^^] = (2p- l)i?[£^+„eP,]. 

The result is then obtained by recurrence. □ 
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We now introduce an extra randomness in the persistence. We first denote by 
PP the law of for a given p. Now, considering a probabihty measure /i on [0, 1], 
we will call , the annealed law of the correlated walk associated to /i, i.e. the 
measure on defined by P^" := j}pPdnip). 

Remark: unlike the situation in Q, the persistence does not depend on the level. 
Only one coin toss, according to /i, decides for the whole environment. 

Let be a process of law P^. Let us now introduce the notation := 
— X^_-j^. From Proposition 1 we get the straightforward result: 

Proposition 2. Vm > l,n > 0, E[e!^e'^+„] = Jq{2p- 

2.2. Statement and proof of the result. The goal now is to introduce a prob- 
ability measure fi leading to the same equivalent as (1), mentioned in the intro- 
duction, so that by taking the average over a large number of trajectories, we 
approximate a discrete Gaussian process having the same properties as in |p^ , 
whose scaling limit is the fractional Brownian motion: 

Theorem 1. Let H e]l/2,l[. 

Denote by the probability on [i, 1] with density (1 - H)2^^'^" {1 - pY^"^^ . 
Let [X^ '*)i>i he a sequence of independent processes of law P^ , 

hm C hm ^ 



With CH = y r(3-2g) ■ 

C means convergence in the sense of the finite-dimensional distributions, and 
means convergence in the sense of the weak convergence in the Skorohod topology 
on -D[0, 1], the space of cadlag functions on [0,1]. 

Proof: The central limit theorem implies that £limjv,/^oo — 



discrete centered Gaussian process {Y^)k>i, with stationary increments : 
- y" "th E[G«] = 0. E\[G«f\ = 1 and 



Q 



r{n) =(2-2i?)/ v'\l - vf-'^" dv 

i)r(2-2g) 

+ i-2H) 

1 1 H{2H ~1) 



f (n + 3 — 2H) n 



r(3-2iJ) 2 - 2 



So that. 



N N 



--1 j=l 

N-1 



4W0)+EH0) + 2^r(fc)]) 

1=1 k=l 
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(the last step consists simply in two successive comparisons between sums and 
integrals). A direct application of (lemma 5.1) allows to conclude. □ 
We can give also an analog statement for _ff = 1/2: 

Theorem 2. Denote by /i^ the uniform probability on [i, 1]. 

1 1 
Let (X^''*)i>i be a sequence of independent processes of law P^', 

lim C lim t^^' ^s{t) 

AT^oo M^oo 2 ^ N log NVM 

where B is the classical Brownian motion, and ci = -j= . 

2 v2 

Proof: The scheme is the same as in Theorem 1. The difference here is that 
r{n)^2jl{2u-irdu^j^. 

So that' r(0) + EZ'i^HO) + 2 ^^^^ r(fc)] 2iVlogiV. 

We conclude again, applying [|3| (lemma 5.1). □ 
Remark: The order of the limits in both theorems is of big importance: the limit 
in the reverse order would bring 0, as far as for any fixed p, a correlated random 



walk satisfies a central limit theorem with normalization \/ N — o{^/N\ogN) (see 
i)- 

We want now to compare previous theorems with the result of |Q : in , the 
limit theorem consists also in taking the scaling limit of the average over a large 
number (tending to infinity) of i.i.d. copies of processes which are a succession of 
falls and rises. In the case of ||l^, falls and rises are all independent with infinite 
variance. 

In our setting, even if the lengths of falls and rises have finite variance under 
each (they are geometric) , the laws under P'^ of the falls and the rises (which 
are all the same) have infinite variance: indeed, if we denote by L the length of a 
rise, 

P^" {L>n) = /p PP{L > n)dpL"{p) = p''dn"{p) 

= {2-2H) f\"il~p)^-'"dp ^ T{3~2H)^^ 

Jq n^oo ri 

(for i? = 1/2, we get P^* {L > n) - -) 



and we get the same kind of tail as in |14[ . 

The difference lies in the fact that the falls and the rises are not independent: 
indeed, if the first rise is short, it probably means that the environment p is small, 
so that the following fall will be probably short also. More precisely, if Li and L2 
denote the lengths of the first rise (resp. fall) and of its following fall (resp. rise), 

P^" (Pi >n,L2>m) = E^H [PP{Li > n, L2 > m)] = E^h 

r(3-2iJ)- ■ .2-2H 

n— ►00 (^n + rn) 

Li and L2 are therefore not independent. 

Note: Theorems 1 and 2 can be extended to any probability measure with mo- 
ments equivalent to ^2-2Hi^f„) , where L is a slowly varying function. The additional 
arguments can be found in p3| and rely mainly on Karamata's theorem in order 
to replace the naive comparison between sums and integrals we used at the end of 
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the proof. We chose /i^ in our statement for the simple formula it yields for ch 

1 

and because it is the law of 1 — ^^-^ — where U is uniform on [0, 1], which makes 
it easy to simulate. At the end of the article, we will discuss the practical interest 
of other measures. 

3. The case < F < 1/2 

3.1. The alternating correlated random walk. We first remark that the cor- 
related random walks of section 2, cannot provide negative correlations for the 
increments, at least for increments separated by an even time interval. The best 
we can hope is to get an alternate sign for the correlations. In order to get a 
process with always negative correlations (except for variances), we will consider 
the sequence of the sum of two consecutive increments. More precisely, if we con- 
sider, with the notations of section 2, the sequence (efn+i + ^^n+2)n>o^ for ^-ny P 
less than 1/2, we get indeed a sequence of negatively correlated variables and it is 
also possible to exhibit a probability on [0, 1] such that the equivalence relation (1) 
mentioned in the introduction will be satisfied. 

But, in the case < < 1/2, this condition alone does not ensure a scaling 
limit for the scale (which is that time smaller than ^/N). It has to be allied to 
a compensation relation between all the correlations. We refer to |l^ (section 5) 
for the statement of this condition, and we will explicit it further. 

It is the reason why we have to introduce a different kind of walk, we will call 
alternating correlated random walk with persistence p. It is a process evolving on 
Z by jumps of -1-1 or -1, whose probability of making the same jump as the previous 
one is alternately p and 0. In other words, one jump over two is the opposite of the 
previous one: 

Definition 2. For any p S [0,1], the alternating correlated random walk with 
persistence p, is a Z-valued discrete process, such that : 
_ XP = 0, PiXf^ -IH 1/2, PiXf = 1) - 1/2. 
_ Vn > 1, £P := XP — Xfi^_i equals 1 or -1 a.s. 
_ Vn > 1, P{el, = el^^MXl < k < 2n - I)) = p. 
Vn > l,eL+i = -^2n- 



As suggested in the introduction of this section, we will be actually interested 
by the process 

v-p 

yp ^2n 

" ■ 2^ 

(the importance of this normalization will appear later). The trajectories of this 
process take only two values, which are either — 1/^y/p and or and Xj ^fp. The 
successive lengths of the time intervals during which the process stays on each value 
are independent and geometric with parameter \ — p. 

We compute now the correlations of the increments of i.e. of (5^ :— — 

Proposition 3. _ Vm > 1, E[{5^f] = 1 

_ Vm > 1,71 > 1, E[6P^5l^^] = -p{l - 2p)"-i 

Proof: As in Proposition 1, everything is based on the following facts: 
Vn > 1, 
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_ E[eP,Ja{Xl, 0<k<2n-l)] = i2p- l)eP,„_„ 
_ i;[£t„+ik(X^,0 < k < 2n)] = 

Vm > 1, 

= ^(2 + 2(2p-l)) 

= 1 

Vm > l,n > 1, by successive conditionings, 

-^['^ni'^TO+n] = 4p-^[(^2m-l + ^m)(^m+2n-l + ^2m+2n)] 

= |(1 - 2p)"-ii?[(eX„_i + eiJ(eL+i + el 
= X(i _ 2p)"-i(-l - 2(2p - 1) - (2p - 1)2) 
= -p{l - 2p)"-i 



□ 



We can already note the following (compensation) relation: 
«n+2E^[«+J=0 (2) 

n>l 

Again we introduce an extra randomness in the persistence. We first denote by 
QP the law of for a given p. Now, considering a probability measure on [0, 1], 
we will call Q^, the annealed law of the correlated walk associated to fi, i.e. the 
measure on defined by dQ'' := QPdij{p). 

Let y be a process of law Q^^. We introduce := ^ — Prom Proposition 
3 we get the straightforward result: 

Proposition 4. _ Vm > 1, ^[(^^J,)^] = 1. 

_ Vm > 1, n > 1, J = - /o P(l - 2p)»-idM(p)- 

Note: As the compensation relation (2) is satisfied "p by p" , it remains true for 
the annealed correlations. 

The goal now is to introduce a probability measure leading to the same equiv- 
alent as (1), mentioned in the introduction. 

3.2. Statement and proof of the result. We proceed as in previous section: 

Theorem 3. Let H e]0, l/2[. 

Denote by the probability on [0, i] with density (1 — 2H)2^~'^^p~'^^ . 
Let (F" '%>i be a sequence of independent processes of law , 

hm £ lim + 



withcH = y r(i-2H) ■ 

Proof: The central limit theorem implies that CVoxim ^oo — ^JW^ ^ 

discrete centered Gaussian process {Z^)k>i, with stationary increments := 
Z^+i - Zk with E[W^] = 0, E[{W^f] = 1 and 



1 

yi,n> 0, r(n) := E[W^"W^^^] = -(1 - 2H)2^-'^" (1 - 2w)V" 

Jo 



2H 



du 
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{l-2H)Tln + l)T{2-2H) 
^ 2 r(n + 3 - 2H) 

W N 



So that, 



i?[cl,(Gf + ... + G^)2] =.4^^r(|z-j|) 

1=1 j=i 



= 4(r(0)+^[r(0) + 2^r(fc)]) 



i=l k=l 
N-1 oo 



= 4(r(0)-2^ 5] r(fc)) 



n — ^oo 



The last equahty comes from the compensation relation (2), and the last step 
consists simply in two successive comparisons between sums and integrals. A direct 
application of (lemma 5.1) allows to conclude. □ 

The process reminds strikingly the On/Off processes of the traffic modeling 



theory, described in or the "1/" -renewal process of 12 which would have 

alternate rewards. But our attempt to make this model work here, failed probably 
because it is not clear how to translate to this model the normalization of by 
Y^, which is crucial in our construction because its absence would yield an infinite 
measure for /i^. 

Remark 1: As in previous section, /i^ is not the only convenient measure. We 

1 

remark that it is the law of ^ ^ ' where U is uniform on [0, 1]. 

Remark 2: Other normalizations of by powers of p can be considered, but 
then the measure fi^ has to be multiplied by a power of p also. 



4. Practical aspects 

As mentioned above, taking the limit in reverse order, yields trivial processes. 
The question is: for a given number of steps N , what kind of M are we supposed 
to take in order to approximate the right limit ? We restrict our study to the law 
of Bh{^)- We base our study on Berry-Esseen's inequality, applied to the sequence 
of i.i.d. variables with law crX'^" /N" (resp. crY^" /N") for H > 1/2 (rcsp. 
H < 1/2). 

What remains to compute, is an upper bound for the third moment of the ab- 
solute value of these variables. 
We begin with the case H > ^. 

Proposition 5. For H > 1/2, for N large enough, 
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Nota Bene: we express Dh in terms of ch, in order to have a formula that works 
for other measures, as it wiU be useful in the following. 

Proof: We omit here the superscipts /i^ in the variables Xj:^ 's and eJJ 's. = 

Y.k=i^k, with: 

_ £„ are Bernoulli(l/2), 

_ cjfCoviek, ei) = r{\k - l\), with r(n) --oo -^^i^s^- 
We use Cauchy-Schwarz inequality to get: 

^2^ 12 

l<ik<N 
l<fe<4 



(using that the variance converges to 1 for large N). 

Assume 14 < «3 < «2 ^ *ij we get, as in Proposition 1, by successive condition- 
ings: 

So that, 

E[ei^ei^Si^ei^] = r{{ii - 12) + («3 - u)) 
We see by using the equivalence between sums and integrals, as we did at the end 
of the proof of Theorem 1, that the sum is equivalent to: 

i-N i-xi /•X2 rxs 

4! X / / / / r{{xi — X2) + (x3 — X4))dx4dx3dx2dxi 



^1 pN rXi i'X2 PX3 

~ — X H{2H - 1) / / / / ((a;i - X2) + {x^ - X4)f'"~'^dx4dxzdx2dxi 

Jo Jo Jo Jo 

_ 6 {2H-l)N^"+^ 
~ 7^ {H + 1)(2H + 1) 

□ 

Now applying Berry-Esseen's inequality, the error of the distribution function of 
the marginal at time 1, is dominated by 0.65 x Dh ^^jj (using the constant 0.65 of 

(Tsl , as far as the third moment is much bigger than the power | of the variance). 

We deduce that Theorem 1 remains true as soon as M{N) ^ 00 as iV ^ 00 at 
a faster rate than N^~^^: 

Corollary 1. Let M he a function on the integers such that ^^a-^l tends to 00, 

£ lim CH — , = OHit) 

JV^oo N"^/M{N) ^' 

Proof: Using the generalization of Berry Esseen's inequality to multidimensional 
variables we get the convergence for finite-dimensional marginals. 

To get the weak convergence, we cannot use directly jl^ as in previous section, 
and we prove the tightness of the family of processes, by checking Billingsley's 
criteria Q (Theorem 15.6): 

Denote by 5^.(i) :^ ^^f^' ■ 
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Let 1 > <2 > t > > 0, and A: e N, 

When we develope the polynomial with degree 2k inside the expectation, we 
notice that only the monomials of the type (X^)^""^ ...(X^)^"*^ have a non null 
contribution, as far as the X*'s are centered and independent. This set has cardinal 
0(M'=). 

Now, we find, by the same means than in Proposition 1, that 
Vai > 1, E[Xl'''] = 0(Ar2"i+(2H-2))^ 

This leads to E[SNit)^''] = 0{{Ntf''+^^"-^^'') = 0{{Ntf"''). 
Hence, for some positive constant C, 
JN{kM,tM) < C{h - tY"\t - hf"^ < C{t2 - 

We choose k > in order to satisfy Billingsley's criteria. □ 
We conjecture that this result remains true in the case of which proposes an 
answer to the question asked at the end of jl^ . 
We do the same for the case H — 

Proposition 6. For N large enough, 

i?[(4£E)3] < c. X 

V^logTV^ ' - 3 logiV 

Proof: As in previous proposition, Xjy = X^fcLi ^k, with: 
_ £„ are Bcrnoulli(l/2), 

_ ciCov{ek,ei) = r{\k - l\), with r{n) --oo ^■ 

Similarly as in previous proposition, ^[( " ) ]^ is equivalent to: 

V N log N 

2 

Ci /J^l pN r-xi PX2 rX3 



But, 



/o Jo Jo Jo i^l - ^2) + {X3 - Xi) 6 

□ 

Now applying Berry-Esseen's inequality, the error of the distribution function of 
the marginal at time 1, is dominated by 1.3 x — ^ — . 

' ''log NVM 

Corollary 2. Let M be a Junction on the integers such that jviog(jv)-^ tends to 00, 



"'0 



NlogN^'c^ J III riixi - X2) + {x3 - X4))dx4dx3dx2dxi)'' 



N nxi /•X2 px^ 

da; 4 da; 3 da; 2 da; 1 = 



2 ^M{N) X N\og{N) 
Finally, we treat the case H < 
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Proposition 7. For H < 1/2, for N large enough, 



with Dh = ^J^Bri ^ 



Proof: The proof is more delicate than in Proposition 5, because the computation 

of the fourth moment makes appear compensations (as in the variance computa- 
tion), and wc have to treat them first "p by p" , before integrating against fj.^ . 

Assume 14 < is < 12 < ii, we get, as in Proposition 1, by successive condition- 
ings: 

where rp(0) = 1 and Vn > 0, rp{n) = — p(l — 2p)"~^. 
Now, 

E =12 E E[61S161SI] + 0{N) 

l<ii^<N max{i3,i4}<min{zi ,12} 

l<fc<4 

= 12 Yl rp{\ii-i2\)rp{\i3-i4\) + 0{N) 

max{i3,i4}<min{n ,12} 
A'' / min{n,i2} \ 

= 12 E E r,{\is-i4\)\rp{\i^-i2\) + 0{N) 

min{ii ,22 } — 1 ymax{i3 ,24 } — 1 / 
TV / min{ii,i2} / max{z3,i4} \ \ 

= 12 E E kp(0) + 2 E rp{l)\\r,{\i,-i2\) + 0{N) 

min{zi ,12}=! \ max{z3 ,24} — 1 \ /=1 / / 

N I min{2i,i2} \ 

= E E (l-2pr""{'^''^> rj,(|ii-i2|) + OW 

min{zi,i2} = l \ max{i3 ,24 }=1 / 
A'^ / min{2i,Z2} \ / N—ui'm{ii,i2} \ 

= 12 E (l-2p)M r,(0) + 2 Y rpil)\+OiN) 

miii{ii,i2}=l V k=l J \ 1=1 J 

N /min{n,i2} \ 

= 51 E (l-2p)M (l-2p)^-™{'i'*=> + 0(^) 

minj^i ,^2 } — 1 \ k=l J 

= 12 Y (1 - 2p)^-('-'=' 0(7V) 

l<k<l<N 

(The contribution of the exceptional situations max{?3,Z4} = min{ii,?'2} is esti- 
mated by 0{N), because of the same compensations as the ones described in the 
above equalities.) 
We now use, 

J^{l-2prdn"{p) =-2^r(fc) 



1 2H 



We deduce. 
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-H 



l<ik<N 
l<fe<4 

(2H /-^ dxi 



N^" \cl Jo Jo {N-{x2-xi)) 



1-2H 



in 1 TT 

□ 

Applying Berry-Esseen's inequality, the error of the distribution function of the 
marginal at time 1, is dominated by 0.65 x Dh^^ 



Corollary 3. Let M he a function on the integers such that ^\^IVh tends to oo, 

N^^ NH^/M{N) ^ ' 

As a general conclusion, we can say that, using /i// , the number of computations 
we have to make, in order to get a N steps trajectory, is a constant times N^~'^^ 
for i < i/, and N^^^^ for < iJ < i. In any case, it is a power of N between 1 
and 2. 

The constant in factor, is of big importance for practical simulation. We took the 
best constant in the Berry-Esseen's inequality we could find in the literature, even if 
it is bigger than the best constant possible conjectured by Esseen, i.e. '^^J^ — 0.41 
that would gain in time a squared factor equal to 2.5 (see [Tsl for a nice discussion 
on this subject). 

But more important is to note that the constant can be considerably ameliorated 
by using other measures than jin, providing smaller ch's. It is the case for the real- 
indexed sequence of probabilities {^iH.k)k>o, defined below: 

-Vffe|i,l]. 

which is just the law of ^+-^(^^^ ^-^) ^j^fj coincides with iih for k = 1 (where B(a, b) 
denotes the Beta variable with parameters a and b). 

_¥Hel0.i|, 

which is just the law of — — 2 — ^^'^ coincides with /iff for fc = 1. 

_ For > 2, we obtain CH,k - ^ r(fc+2-2g) ~ 1^ — ■ 
So that, using ^H,k for large k, yields an error equivalent to: 



^ {H + l)(2H+l) 'k' Vm 
For H — i, we obtain a i. = -^r- 
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The error is equivalent to: 



0.65 ,N,, 



logAT ^fc^ VM 

r,fc 

The error is equivalent to: 



For < i, we obtain CH,k = ^ r(l+l-2H) ~ 



Using ^H,k instead of hh, allows to gain for M a factor (resp. k^~'^^) 

for i? > i (resp. for < ^). Loosely speaking, it erases the noise generated 
by the parameters between | and a fixed constant smaller than 1. The trouble 
making k increase, is that it damages the value of the covariance of Xj^ (resp. 
Yat), but wc can allow any k{N) — o{N). So that we obtain an algorithm with 
any M{N) = l/o(l) number of trajectories. As a result, our algorithm requires a 
number of computations of the order N/o{l), for any o(l). Moreover, we have only 
M{N) real datas to keep in memory along the whole procedure, corresponding to 
the coin-tossed parameters of the walks, and M{N) integers (-1 or 1), giving the 
last moves of the walks. 

We want to give now a second family of measures (/x^ k)k>o- 



_ For if > i, is the law of 1 - i^-u^)™ _ 

An easy computation gives = (error: 0.65 x Dh ^^j^^ 



_ For H < ^, n'jj ^ is the law of (LJiHZZEl. 

Again, an easy computation gives c'^j ^. = ^ (error: 0.65 x Dh ^J^j^ )■ 

The advantage of this family is obviously the easy simulation it provides. The 
error is estimated by a term containing that seems to be better than the last 
one, but the damages on the variances grow faster than for iiH,k- Actually the scale 
\/fc corresponds, in the previous family, to the scale k^^^ (resp. k^^^) for H > j 
(resp. H < i). The drawback of this family, is that the theoretical computations 
of the variance are not very explicit. 

We illustrate our results by three graphs corresponding to three different param- 
eters of H, with N = 1000, and a theoretical error smaller than 10% (we indicate 
the time it takes for Matlab to draw a graph): 

_ For if = 0.25, wc take M = 200 and fc = 1 (15 seconds) (Fig. 1). 

_ For H = 0.5, we can see that the convergence in Theorem 2 is the slowest one, 
and we will use the simple random walk to simulate it. (Fig. 2) 

_ For = 0.75, we take M = 400 and k = 0.5 (25 seconds) (Fig. 3). 

Remark: The interest of taking large k appears when N is very large, especially 
if we want to preserve a small error on the whole trajectory. In the case H > ^, 
we were even obliged to take k smaller than one. 

We remark the different behaviours of the trajectories: let us recall that the 
Hausdorff dimension of the trajectories are a.s. equal to 2 — iJ, and we notice that 
the variances of the process between and 1 become larger when H decreases. 

As we noticed just above, we may be limited by the fidelity of the covariance of 
our process. In this spirit, it is quite interesting to note that the auto covariance 
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function of the Gaussian noise Gnij) of the introduction is, up to a shift, the 
sequence of moments of a probabihty measure on [0,1]. I first remarked it by 
checking the conditions of Hausdorff theorem p. 226), but Marc Yor gave me 
kindly, the method to get the exphcit density of this measure, and I present it here. 
This brings a third (the last !) family of probability measures: 

Proposition 8. Let H 1[. Consider the family of probability measures {vH,k)k>n 
on with density C{H,k) x (1 - p)2(2p - l)'=-i(ln(2^))-i-2H^ 

C{H, k) := iff^l^ X {{k + 2f" - 2{k + + A;2«V' 

{n + k + 2f" - 2(n + fc + 1)2^ + {n + kf" 



Vn>0, / {2p-lTdvH,k{p) = 
Jo 



(fc + 2)2H - 2(A: + l)2-f^ + 
Proof: Vn > 1, 

+ i)2ff _ 2n^" + in- if") [\n + tf""-^ ~ (n + t - lf"-^dt 

2 Jo 

r-l /-l 



= H{2H-l) [ [ {n + t + s - l)^"-^dsdt 
Jo Jo 

T{2-2H)jJJ, 

H{2H - 1) /•+°° ^_„„^^ _ g-«)2g«^-i-2H^^ 



r(2 - 2H) Jo 

^Hi2H-l) r,»(l^).(i,l)-i-2H,^ 
r(2-2i7) Jo X ' ^ x' 

The end of the proof is straightforward, by change of variable. □ 
Note: The normalization corresponding to i'H,k is: 

_ f ik+2)^"-2{k+lf"+k^" ^1 VH{2H-1) 

The same can be done for H e]0, |[: 



Proposition 9. Let H s]0, i[. Consider the family of probability measures [v^ k)k>o 

'2pJ 

V 1(1. _^U2H _ i,2H\-l 

r(l-2H) 

{n + k + 1)2^ - 2{n + kf" + {n + k- if" 



on [0, \], with density C{H, k) x p{2p - if^^iMj^))"^^'^", and 
C{H, k) := X {{k + If" fc2H)-i 



yn>l,-J^pi2p-ir-^d.,Ap)- 2((fc + l)2--P-) 

Proof: Vn > 1, 



i((n + If" - 2n^" + (n - l)^^) =H I (n + tf"~' - {n + t - lf"~'dt 
H 



T{l-2H)Jo Jo 
H 



1 />-\-OQ 

(e-("+*)" - e-^^+'-^^^)u-^"dudt 



r(i - 2H) 



H f ^n(L^)2(i^l)-l-2ff^^ 



r(l - 2i7) Jo ' X ' ' x' 
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We find the normalization, using + fc + 1)^-^ — 2(n + fc)^^ + (n + Zc — 1)^^ = 

k"^" -{k + 1)^", and the relation 1 = 2^„>ip(l - 2p)"-i. □ 
Note: The normalization corresponding to i^H,k is: 

These measures i'H,k have an interest, as far as X'^ (resp. Y^^") have explicit 
covariances matching quite well with the covariances of the fractional Brownian 
motion. On the other hand, they don't seem very easy to simulate. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 1: H=0.25 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 2: H=0.5 




Fig. 3: H=0.75 
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